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Abstract 

The Nordstrom- Vlasov system provides an interesting relativistic generaliza- 
tion of the Vlasov-Poisson system in the gravitational case, even though there is 
no direct physical application. The study of this model will probably lead to a 
better mathematical understanding of the class of non-linear systems consisting 
of hyperbolic and transport equations. In this paper it is shown that solutions of 
the Nordstrom- Vlasov system converge to solutions of the Vlasov-Poisson system 
in a pointwise sense as the speed of light tends to infinity, providing a further and 
rigorous justification of this model as a genuine relativistic generalization of the 
Vlasov-Poisson system. 



1 Introduction 

Kinetic models of collisionless matter have many important physical applications. In 
astrophysics, for example, the stars of a galaxy are often modelled as a large ensemble 
of particles in which collisions are sufficiently rare to be neglected. The distribution 
foo °f particles in the phase-space satisfies the Vlasov-Poisson system: 

fi t / O o+p-V B / oo -V B C/'-Vp/oo=0 I (1.1) 

A X U = An-fpoo, 7 = 1, poo = / foodp. (1.2) 

In the previous equations, foo = foo{t,x,p) gives the probability density to find a 
particle (star) at time t at position x with momentum p, where t € R, x £ R 3 , p € R 3 . 
U = U(t, x) is the mean Newtonian potential generated by the stars. 

*E-Mail: mg026@math.chalmers.se 
^E-mail: hayoung@aei.mpg.de 



1 



THE NORDSTROM-VLASOV SYSTEM 



2 



By replacing 7 = — 1 in i|1.2[l one obtains the Vlasov-Poisson system in the plasma 
physics case. Here the particles are charges and U is the electrostatic potential which 
they create collectively. We consider a single species of particle in both cases. The 
applications of these Vlasov-Poisson systems are restricted to the situations where the 
relativistic effects are negligible, i.e., low velocities and weak fields. Otherwise the 
dynamics has to be described by the relativistic Vlasov-Maxwell system in plasma 
physics and by the Einstein- Vlasov system in stellar dynamics. 

The two Vlasov-Poisson models are very similar to each other and no substantial 
difference arises in the question of global existence of classical solutions, which is by 
now well-understood (cf. O EB 1201 EE!)- As opposed to this, the relativistic models 
have very different structure and so far they have been considered separately. In the 
gravitational case, global existence of (asymptotically flat) solutions for the Einstein- 
Vlasov system is known only for small data with spherical symmetry For the 

relativistic Vlasov-Maxwell system the theory is more developed, cf. |H10> |5]-|15|. 
\'21\ . However global existence and uniqueness of classical solutions for large data in 
three dimensions is still open. 

In a recent paper 0, a different relativistic generalization to the Vlasov-Poisson 
system in the stellar dynamics case has been considered, in which the Vlasov dynamics 
is coupled to a relativistic scalar theory of gravity which goes back, essentially, to 
Nordstrom ^2]. More precisely, the gravitational theory considered in |2] corresponds 
to a reformulation of Nordstrom's theory due to Einstein and Fokker (see 8 ). The 
resulting system has been called Nordstrom- Vlasov system and reads 

- dU + c 2 A x cb = 4tt / fdp (1.3) 
J v 1 + c p 

d t f+p-V x f-[S(cftp+-fid=]-V p f = AS{<i))f, (1.4) 

where 

v 2 = \v\\ p = (i + c- 2 P 2 r 1/2 p, s = d t +p-v x . 

Here / = f(t,x,p), <j> — <f>(t,x) and c denotes the speed of light. A solution (/, <j>) of 
this system is interpreted as follows. The spacetime is a Lorentzian manifold with a 
conformally flat metric which, in the coordinates (ct, x), takes the form 

9llv = e 2 *diag(-l, 1,1,1). 

Throughout the paper Greek indices fj,, v and a run from to 3 and Latin indices a 
and b take values 1, 2, 3. The particle distribution / defined on the mass shell in this 
metric is given by 

f(t,x,p) = e- 4 *f(t,x,e*p). 

More details on the derivation of this system are given in the next section. It should 
be emphasized that, although this model has no direct physical applications, scalar 
fields play a major role in modern theories of classical and quantum gravity. For ex- 
ample, the Brans-Dicke gravitational theory yQ, which is continuously tested against 
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general relativity, is a combination of Einstein's and Nordstrom's theory. The Nord- 
strom- Vlasov system is also interesting in a pure mathematical sense. A hope is that 
by studying this model one may reach a better understanding of a class of systems 
consisting of hyperbolic and transport equations. 

However in order to justify this model as a genuine relativistic generalization of the 
(gravitational) Vlasov-Poisson system, it is necessary to indicate the relation between 
the solutions of the two systems. The main goal of this paper is to prove that in the 
non- relativistic limit c — > oo the solutions of 1|1.3[1 - (|1.4|) converge to solutions of Ql.l|) - 
<|1.2[l in a pointwise sense. The analogous result was proved in for the relativistic 
Vlasov-Maxwell system (see JB] for the case of two space dimensions) and in |2Sj for 
the Einstein- Vlasov system with spherical symmetry (in the latter case a weaker form 
of convergence holds also in the absence of symmetries, see |24j ) . 

This paper proceeds as follows. In section El we provide a formal derivation of the 
Nordstrom- Vlasov system and state our main results in full detail. The first of such 
results is a local existence theorem of solutions of the Nordstrom- Vlasov system in an 
interval of time independent of the speed of light, which is a necessary step to proceed 
further in the study of the non-relativistic limit. The solution of the latter problem is 
our second result. Our analysis follows to a large extent and is based on the use 
of certain representation formulae for the solutions of the Nordstrom- Vlasov system 
which have been introduced in [3] and which will be adapted to the present case in 
section [3] There we shall also prove some estimates needed in the sequel. One of 
these estimates states that the distribution function / is uniformly bounded, which 
permits to improve the conditional global existence result, Theorem 1 in 0]. (The 
general L q estimates on / are proved in [2], where they are used to establish existence 
of global weak solutions to the Nordstrom- Vlasov system). In section 4 we prove our 
main results. 



2 Derivation of the Nordstrom-Vlasov system and 
main results 

We shall refer to the Nordstrom-Vlasov system as the set of equations which models 
the kinetic motion of a self-gravitating ensemble of collisionlcss particles in accordance 
to a gravitational theory satisfying the following assumptions: 

(1) The gravitational forces are mediated by a scalar field (j> and the effect of such forces 
is to conformally rescale the metric of the (four dimensional) spacetime according 
to the relation 

g = A 2 (<p) v , (2.1) 
where n is the Minkowski metric and A is a positive function. 

(2) Scale invariance property: There exists a one-parameter symmetry group whose 
action consists in rescaling A((j)) by a constant factor. 

(3) Postulate of simplicity: The dynamics of the field cf> is governed by second order 
differential equations. 
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(4) The matter (by which we mean any non-gravitational field) is universally coupled 
to the metric Q2.1|) . 

It was observed in [2j (appealing to the more general case of Scalar-Tensor theories 
considered in [H]) that the above assumptions single out a unique one-parameter family 
of scalar gravitation theories. This parameter appears because of the scale invariance 
property which forces the conformal factor to be of the form A((j>) = exp(K0), with 
k > 0. Hence in this theory the spacetime is a Lorentzian manifold endowed with the 
metric 

.9 = e 2K V (2.2) 

To write down the field equation of this scalar gravitation theory in a simple form, let 
us consider a system of Cartesian coordinates {a; := ct, x , x 2 , x 3 } — {x^}, c denoting 
the vacuum speed of light in Galilean frames. In these coordinates, the equation for <f> 
takes the form 

-c- 2 d?<j) + A x <f) = -47r% K e 4K<# T. (2.3) 
c 4 

Here G* is a dimensional constant (the bare gravitational constant) and T is the trace 
of the stress-energy tensor of the matter with respect to the physical metric g. In 
12] c, G» and n have been set equal to unity and the factor in has been removed for 
simplicity. As we already mentioned in the introduction, this scalar gravitation theory 
corresponds to the one considered in [HI . 

In the case of the Nordstrom- Vlasov system, the dynamics of the matter is described 
by a non-negative, real-valued function / which gives the probability density to find 
a particle in a given spacetime position x^ and with a given four momentum p^ . We 
assume for simplicity that there is only one species of particle and choose units such 
that the proper mass of each particle is equal to one. The particle distribution / is 
defined on the mass shell of the metric (|2.2() , which is the subset of the tangent bundle 
of spacetime defined by the condition g liv p ti p v = — c 2 , p° > 0. This implies 



p" = Ve- 2 ^ 2 + 6 ab p a p b . (2.4) 

Using (x ll ,p a ) as coordinates on the mass shell and denoting by dp the volume element 
dp 1 dp 2 dp 3 , the stress-energy tensor for this matter model is 

T"" = -cf V\tetg\—fd P , 

J Po 

which implies 

T= -c 3 e 2K<?i [ ^dp. (2.5) 
J P' 

Finally, the coupling between the scalar gravitational field and the matter is completed 
by requiring that the distribution / of particles on the mass-shell is constant on the 
geodesies of the metric (|2.2f) . This leads to the Vlasov equation: 

c-'dj+^f- ^T« u d pa f = 0, (2.6) 
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where T° — n^d^cf) + S^d^cf) — r]^ u d a (f)) are the Christoffel symbols of the metric 

Our goal is to relate the solutions of the system (|2.3[) - H2.6[) to the solutions of the 
Vlasov-Poisson system <|1.1[1 - H1.2ll satisfying the condition limi. I .|_ i . c>0 U — (isolated 
solutions). Hence (/oo,J7) solves the system 

dtfoo + V ■ V./oo - V X U • Vp/oo = 0, (2.7) 
U=-G f^^-dy, (2.8) 

J \y-x\ 



p ao (t,x) = / f oc (t,x,p)dp, (2.9) 

where G denotes the Newtonian gravitational constant which had been set equal to 
unity in (|1.1|) - (|1.2I) . In order to get some light on the relation between the two systems, 
let us consider a formal expansion of the solutions of the Nordstrom- Vlasov system in 
powers of 1/c: 

<j) — 00 + C _1 0i + C~ 2 (f>2 + ... 

.f = h + c-\h+c- 2 h + ... 
Substituting these into l|2.3(l and comparing the terms of the same order we obtain 

A^ = 0, A x 0i = 0, (2.10) 
-d 2 (j} + A x (j) 2 = 47rGUe 7K ^ J f Q dp. (2.11) 

Assuming fields vanishing at infinity, (|2.1U|) implies <fio = <fii = and so l|2.11|l reduces 
to (|2.8|l with the identification (f>2 ~ U, /o ~ foe, provided that G*k = G. The latter 
condition, which is necessary in order to obtain the correct Newtonian limit, shows 
that the role of the scale invariance parameter k is merely the one of fixing the units 
of the corresponding theory. We shall henceforth set n = G* = G = 1 for simplicity. 

To put the above formal discussion in a more rigorous mathematical context, we 
first rewrite the equations H2.3() - (|2.6() with the "unphysical" particle density as in the 
formulation of 4 , namely 

f(t,x,p) = e^f(t 7 x,e-*p). 

In this frame, the unknown (/, cj)) satisfies the equations l|1.3|l - H1.4|) . We supply this 
system with initial data < f(0,x,p) = f n {x,p), 0(0, x) = 4>^{x), d t (p(0,x) = ^"(x). 

The following notation will be used. Given two functions g and h on M™ we write 
g < h if the estimate g < Dh holds for a non-negative constant D independent of 
c > 1. The constant D may also depend on the length of some time interval [0,T], in 
which case we write g < h for t € [0, T]. Furthermore we write 



A = B + 0(c- s ), 5>1, 
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if \A(y) - B(y)\ < c' 3 , Vy £ R". We also set 

V c {t) = sup {\p\ : (x,p) G supp /(«)} + 1, (2.12) 

0<s<t 

where supp f(t) means the support of f(t, x,p) on (x,p) G R 6 for each i. 
Here are the main results of this paper: 

Theorem 1 Initial data / in G C^R 6 ), ^ G C b 3 (R 3 ), 0i n G C 2 (R 3 ) launch a unique 
classical solution (f,<fi) G C 1 ([0, T max ) x R 6 ) x C 2 ([0, T max ) x R 3 ) fo the Cauchy problem 
for the Nordstrom- Vlasov system IJ..^) - f77^| ) m a maximal interval of time [0, T max ). 
IfPc{Tmax) < °Ot i/ien T max = oo, i.e., the solution is global. 

Note that under the assumptions of Theorem ^ the local time of existence may 
shrink to zero as the speed of light tends to +oo. To remove this possibility we specify 
more restrictive initial data: 

Theorem 2 Assume f ia G C 3 (R 6 ) and <p™ = c~ 2 g^ <£f = C - 2 h\ where 

gKx) = -JJ^^-d P dy, ft«£C c 2 (l 3 ). 

Corresponding to these data there exists a unique solution (f,<fi) € C 1 ([0,T) x R 6 ) x 
C 2 ([0, T) x R 3 ) of f-Z.,'^) - f7?7| ) in an interval of time [0,T) independent of c such that 
Vc(T)<l. 

Finally we give the conditions under which solutions to the Nordstrom- Vlasov sys- 
tem converge in a pointwise sense to solutions of the Vlasov-Poisson system in the 
non-relativistic limit. 

Theorem 3 Let the data for the Nordstrom- Vlasov system be given as in Theorem \^ 
and assume that 

(*) There exists a unique solution (f,<j>) G C 1 ([0 ) T) x R 6 ) x C 2 ([0,T) x R 3 ) of 
\ 1.3\) - \l1\i in an interval [0,T) C [0, +oo) independent of c and the estimate 
Vc{T') < 1 holds for all < V < T. 

Denote by /oo G C 1 ([0, oo) x R 6 ) the global solution of p!.7| )- ifl!.£j) with data f in , which 
is known to exist by ]iy( . Then for every T' G [0, T) and t G [0, T']: 

d t (j){t) = 0{c- 1 ), c 2 ^t) = U(t) + 0(c- 1 ), c 2 V x <t>(t)=V x U(t)+0(c- 1 ), (2.13) 

/(<) = /co(<) + 0(c- 1 ). (2.14) 

In the notation of the spaces of functions used above, the subscript c indicates 
that functions are compactly supported and b means that all the derivatives up to the 
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indicated order are bounded. Theorem^will be proved as a corollary of Theorem 1 of 
£U in the next section and Theorems [21 and [21 will be proved in Section 01 

To conclude this section, we remark that more general data for the field are allowed 
in Theorems [21 and [31 For instance we may require that (f>™ = c~ 2 g$ + c~ 3 g b , where 
g € (M 3 ). This modifies our estimates only by terms which are of higher order in 
powers of c , without affecting the general argument. In order to make the estimates 
below more transparent, we take for simplicity g b = 0. 



3 Preliminaries and proof of Theorem [T] 

First note that the classical solution of (|1.3fl is 

f(t - c- 1 ^ - x\,y,p) 



(f>(t, X) = </>hom(i, X) j 



y — x\<ct 



yjl + c 2 p 2 \y - x\ 



dpdy 



■= 4>hom(t,x) +lp(t,x), 



(3.1) 



where 



cphom (*, x)=d t {^- J 4" ( x + ctuj) dw) + ^- J <f)\ n (x + ctuo) 



du (3.2) 



is the solution of the homogeneous wave equation with data </>Q n and cj)™ and ip the 
solution of l|1.3|) with trivial data. We start with estimating the homogeneous part of 
the field tj>. 

Lemma 1 Let the initial data for the field be given as in Theorems^ and\^ Then we 
have 

||<?(>hom(t)||oo <C-\l+t). 

Proof: By means of (|3.2[) and the assumptions on the data <f>™ and 4>f- we have 



</>hom(c H,X) 



t 

V47T 



g\x + tuS) dui 



c-H 

47T 



fv(x + tu>) du 



(3.3) 



The term in the square brackets in (|3.3I) is estimated by D(l + t) where 

-D := Halloo + HV^Iloo + ll^lloo- 

Hence ||0ham(c _: 4)||oo < cr 2 (l + t), which implies ||0hom(*)||oo ^ c _1 (l + t). □ 

The following estimate is crucial for extending the argument of [35] to the Nord- 
strom- Vlasov system. 

Proposition 1 The distribution function f satisfies the estimate 
ll/WIU < ll/ in ||oo exp [4(110^(011^ + HCIU 
In particular, for data as in Theorems^ and\^ we have ||/(t)||oo ^ e iDt for all t € 
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Proof : Let (X, P)(s,t,x,p) denote the characteristics of (|1.4f) which satisfy the con- 
dition (X, P)(t,t, x,p) — (x,p). In short, we use X(s) := X(s,t,x,p) and P(s) := 
P(s,t,x,p) instead. Note that the function e~ 4 ^f is constant along these curves. 
Hence the solution of l|1.4(l is given by 

f(t,x,p) = / in (X(0),P(0))exp [4<t>(t,x) -40 in (X(O))] 

= / in (X(0), P(0)) exp [-4<£ in (X (0))] exp [4</> hom (i, a:)] exp [4^(i, x)] . (3.4) 

Since ip < 0, then e 4 ^ < 1 and the claim follows. □ 

Combining this result with the one in 0] we obtain the following. 
Proof of Theorem^; It is enough to prove the theorem for c = 1. For given a so- 
lution (/, (f>) of Ql.'ljl - qi.4|l . then the solution obtained by the rescaling c/(c _1 i, x, cp), 
(f>(c~ 1 t,x) solves the system with c = 1. The claim has been proved in £Q under the 
additional condition that Q{T mayi ) < oo, where 

Q(t) = sup {\<f>(t,x)\ : (x,p) S supp/(s)}. 

0<s<t 

Now, assuming compact support in p for the distribution function, i.e., 7 : 'i(T max ) < oo, 
it follows by Proposition that the right hand side of H3.1JI is bounded in L°°. Hence 
(j) itself is bounded and thus the condition <2(T max ) < oo is satisfied. The claim follows 
by Theorem 1 of HJ. □ 

Next we derive the representation formulae for the first order derivatives of the field 
as in 0], but for arbitrary values of c and also with all data terms specified. 
One can see by the Vlasov equation l|1.4(l that 

Sf = {{d t 4>)p + c 2 (l + c-V r 1/2 V^} . V p / + i(S<P)f. 

From H'6. lfl we have 

d t (p(t,x) = <9 t 0hom(i, x) - c -2 ^ 1 / f t ( y,P ^=dpdS y 



v \y-x\=ct J y/l + C V 

- c- 2 [ [ d t f{t-c^\ y -xl m p) dpdy _ 

J\y-x\<ctJ \/l + C 2 p 2 \y — x\ 

Now using the identity 

d t g(t - c-^y - x\,y,p) - (1 + c~ x u ■ p)- x {{Sg)(t - C - l \y - x\,y,p) 

-P- ^y[9{t - cT x \y - x\,y,p)]} 

and integration by parts we achieve the following representation for dt4>- 
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Proposition 2 

d t (f>(t,x) = d t <t>hom(t,x) 

c -w I 1^4 dPdSy 

J\y-x\=ctJ (1 + C X U3 ■p)y/l + C 2 p 2 



2 I I a**(uj,p)f(t- c-^y - x\,y,p)dp-^-2 

J\v-x\<ctJ \y-x\ 



C 

'\y-x\<ct 



f f b**(u;,p)(S<j>)f(t-c- 1 \y-x\,y,p)dp-^— 

J\y-x\<ctJ \V - x \ 

-c- 1 f f c*<(u;,p)(V x 4>)f(t-c- 1 \y-x\,y,p)dp- °" 

J \y-x\<ct J 



i\ y -x\<ctJ \y-x\ 

where the kernels are 



p-(ui + c 1 p) 


(1+ 


C 1 L0 ■ p) 2 \Jl + C 2 p 2 




(lo + c^p) 2 


(1 + C" 


~ 1 ui ■ p) 2 \/l + c~ 2 p 2 ' 




LO + C~ X p 


(1 + C" 


"^•^(l + c-V) 372 



b+ t (u>,p) = 
c*>(u,p) = 
u = {y-x)/\y-x\. 

The process to obtain the representation for V x <t> is similar to the way for d t <f>, but 
now with the following identity: 

d Vi g{t- c _1 |y - x\,y,p) = c _1 Wj(l + c _1 w • p) -1 ^)^ - c _1 |y - x\,y,p) 

+ ( S '' k ~ i + c -i^.fi ) dyk [ 9 ( t ~ c ~ 1 \ y ~ x \' y ' p ^' 1 = 1 ' 2 ' 3 - 

Proposition 3 The representation of d Xi <fi is 
d Xi <j)(t,x) = d Xi (f> hom (t,x) 

-c-H- 1 [ [ ■.f n (y,p)dpdS y 

J\y-x\=ctJ (1 + C 1 W -p)y/l+C 2 p 2 

-c- 3 [ [ a^(cj,p)f(t-c- 1 \y-x\,y,p)dp dy 

J\y-x\<ctJ W X \ 

-c- 3 f f b^(u,p)(Scl>)f(t-c- 1 \y-x\,y,p)dp r ^ 

J\y-x\<ctJ W ~ 



X 



-c- 2 f f c^{uo 1 p){V x ^.f{t-c- 1 \y-x\,y 1 p)dp- 

J\y-x\<ctJ \y-X\ 



dy 
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where the kernels are 

c(uj + c~ 1 p) l - c- l {p A {ujAp))i 



o*-i(«,p)= . 

(l + c- 1 UJ-p) 2 y / l + C- 2 p 2 

C^i (ti),p) = LOiC^. 

4 Proof of Theorems [2] and [3] 

In this section we prove our main results. We shall frequently use Lemmas 1 and 2 of 
|25| . which we state below for future reference. 

Lemma 2 For all g G C^(IR 3 ), we have 



t \g(x + £u)\du<l, 

J\u\=l 

for £ > 0. 

Lemma 3 Let h G C 2 (R 3 ) suc/i t/iaf A/i G C C °(K 3 ). Then for c> and £ > 0, 

A%) 



ft U 



ft / h(x + ctu) dJ) =- -, — dy 
V 7i,.,i_i / iy — sc 



|u)|=l ' — x|>ct I?/ ^1 

The next Lemma contains two simple estimates which are often used in the sequel. 

Lemma 4 (i) (1 + c~ 1 ui ■ p)^ 1 < V c (t) 2 , for (x,p) £ supp fit) 

(ii) f(t,x,p) = 0, for \x\ > R + V c (t)t, where R := sup{|a;| : (x,p) £ supp / ln }. 
Proof: By \p\ < V c {t), 

\ I p>l ,' P ' = , f > 



tJc^Tp 2 ^Tp^i^Tp 2 + \ P \) " 2(c 2 + v c (t) 2 ) ' 

by which (i) follows. The property (ii) on the support of / follows by l|3.4[) and the 
definition of characteristics. □ 

Proof of Theorem^ From Proposition |31 we have 

d Xi (j>{t,x) = d Xt cj) hom (t,x) 

-c-'t- 1 I I g -.r(y,p)dpdS y 

J\y-x\=ctJ (1 + C L LU-p)^l + C £ p z 

+ I Xi +II Xi+ III Xi , (4.1) 
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t 

s47T 



d Xi g Hx + ctu>) du) 



where 

d x ,(/)hom(t,x) = c~ 2 d t (- 

By the assumption £ C 2 (R 3 ), Lemma gives 

c-H 



in 



d Xi h\x + etui) dui 



\w\=l 



in 



<rr 3 . 



d Xi hP(x + etui) du. 

(4.2) 



For the first term in l|4.2|) . using Lemma l^and the identity Ag" = in J f ln dp we get 

< /• ~ ». . . . \ ~ r I f m (v,p) 

l\y-x\>ctJ \y~x\ 



c?tf-; — / d x .g^(x + ctuj)dJ) = -d Xi 

J\u\ = l / 



(ay 



y— x\=ct 



\y-x\>ct 



■ dpdy 



u t f m (y,p)dpdSy 
— 12~ dpdy. 

\y - x\ z 



The surface integral will be combined with the second term in l|4.1|) . First note that, 
using (i) in Lemma 0] and \p\ < V c (t), 



1 



< 



(1 + c- x u -p)\Ji + c~V 

1 



1 



1 



1 



< c- 2 ^ c (i) 2 + 



1 



c _1 |cj • p\ 



c 2 p 2 



yji + C- V 



(1 + c- 1 w -p) 



< c -Xt> 



(4.3) 



So we get 



|j/-x|=ci 



Uif m (y,p)dpdS y 



UiP n {y,p)dpdS y 



y-x\=ctJ (1 + c ^ • p)yl + c 2 p 2 



< c-h- 1 / / C - 1 P c W 3 |^|r(y,p)dp^ y 

«/ \y — x\—ct - 



<c- z V c {tf(ct j f in (x + ctu,p)dpduj' S j <c~ 3 V c (t) 3 



Split the kernel in I Xi according to 
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By (i) in Lemma 0] one can see that \a? 
Lemma 0] we obtain 



,-3 



\y-x\<ct 



a^*i{u),p)f{t- c l \y - x\,y,p)dp 



T^cit) 5 - Using Proposition and (ii) of 
dy 



\y-x\ 



< c 



\v\<R+v c (t)t J\p\<r c (t) 



|a 0Xl (w,p)|/(t- c 1 \y-x\,y,p)dp 



dy 



\y - A 2 



<cr 3 sup Wf^WMt) 5 [ [dp 
o<r<t J \y\<R+tv a (t) J\p\<r c (t) 



dy 



\y - x\ 2 



<c- 3 V c (t) 9 (l + t) 



iDt 



A computation similar to l|4.3|) shows that 

1 



1 



(1 + C^LO ■ p) 2 a/1 + c~ 2 p 2 



< c- l V c {tf 



Hence 



/ / Uif(t~ c x \y- x\,y : p)dp dy 

J\y~x\<ctJ \U ~ x 



\y—x\<ctJ (1 + C 1 UJ ■ p) 2 y/l + C 2 p 2 

f „-l-n W5 

Jh 

<c-' A V c {tf{\ + t)e 
Combining the estimates obtained thus far we get 



f(t-c 1 \y - x\,y,p)dp 



dy 



<c- 2 / / c- 1 V c (t) 5 \u H \f(t-c- 1 \y-x\,y,p)dp 

J \y\<R+V c (t)t J\p\<V c (t) 
3t> /+\9/i , .\ 4Dt 



\y - A 2 

dy 
\y-x\ 2 



-3 + -l 



d Xi 4)hom{t,x) - C t 



y-x\ 



=ctJ (1 + C 1 Ll> ■ p)\/l + C 2 p 



: r(y,p)dpdS y +I XH 



< 



i/(max{0,t-c l \y - x\},y,p)dp 



dy 



dy 



\y-x\< 



+ c~ 6 V c {tfe 



9„4Dt 



^c-^t)" sup H/Wlloo / j^ + c-'Vaitfe 

[0,t] J\y\<R+V c (t)t 

<c- 2 T c (t) 9 (l + t)e 4Dt . 



9 4Z?t 



\y\ 2 



(4.4) 



Now we estimate II Xi and III Xi . Again by (i) of Lemma 0] we obtain the bounds 
\b^(cu,p)\<V c (t)\ | c ^,p)|<7> c (i) 4 - 

Let us define 

K c (t) = sup{c\d t <}>(t,x)\ + c 2 \V x 4>{t,x)\, x S R 3 }. 
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Hence, using c|5(0)| < K c (t), we get 



\H Xi \<c- s V c (t) 4 [ [ f fir.y^y-x^K^dpdSydT 

JO J\y-x\=c{t-r) J\p\<V c (t) 

< c- 2 r c (t) 7 te 4Dt f K c (r)dT. 
Jo 

III Xi satisfies an identical estimate, since c^V^I < K c (t). Collecting the various 
bounds we obtain 

|V„^(t,s)| <c- 2 V c (t) 9 (l + t)e iDt (l + J* K c (T)dry (4.5) 
Now the estimate on the time derivative of the field. From Proposition |3 we have 

f a (y,p) 



d t (j)(t, x) = <9 t </>hom(i, x)-c t 



y-x\=ctJ (1 + C 1 LU ■ p)y/l + C 



2p2 



dpdSy 



I t + lit + 1 lit, 



(4.6) 



where 



<%</>hom(M) 



d t (- — - f h t (x + ctuj) dm) + — [ Ag s (x + ctu) du; 



A direct application of Lemma[21now gives ||9*0hom(*)||oo S c + For the second 
term of (|4.6I) we have 



c-H- 1 



f n (y,p) 



\ y -x\=ctJ (l + c- x tjj ■p)yT+ c- 2 p 2 



dpdSy 



<c-^ c (0) 2 



ct I / f ln (x + ctu/,p)dpdu> 

M = l J\p\<Va(P) 



In order to estimate the remaining terms in l|4.6|l we use the bounds 

\a**(uj,p)\<r c (t) 5 , \b^{u,p)\<V c {tf, \c^(cu,p)\<V c (t)\ 
The estimate for I t follows: 



\It\<c~ 2 V c {tf 



f(t - c 1 \y-x\,y,p)dp 



dy 



\y-x\<ctJ\p\<V a (t) 



|y-x| 2 



< 



c - 2 r c (tf supwmu f i| <c-*v c {tni + t)e iDt . 

[o,t] J\v\<R+v a {t)t \y\ 



Also 



\IIt\<c~ 2 V c (t) 4 / / f{T,y,p)\y-x\-^K c {T)dpdS v dT 

JO J \y-x\=c(t-r) J\p\<V c (r) 

< cr x V c (t) 7 te ADt f K c (r)dT. 
Jo 
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III t satisfies an identical estimate, since c 2 |V x (^>| < K c (t). Collecting the various 
bounds we obtain 



\M(t,x)\ <c~ 1 V c (t) 9 (l + t)e 4Dt (l + j K c {r)d?j. (4.7) 

Combining 14.7|l and 14.5J1 entails 

K c (t) < T c (t) 9 (l + t)e 4Dt (l + J K c (r)dr) . 

Hence by Gronwall's inequality, 

K c {t) < V c (t) 9 (l + t)e iDt exp (V c (t) 9 (l + t) 2 e 4Dt ) . (4.8) 
Note that the characteristics (X,P)(s) of (|1.4|l with (X,P)(t) — (x,p) satisfies 



ds v ™ ' x/TTP 1 

So t t 

\P\<\P(0)\+ [ IC c (T)V c (T)dT <V C (Q)-1+ f IC c (T)V c (T)dT. (4.9) 

Jo Jo 
Therefore by (|4.8|l and the definition of V c (t), (|4.9|) becomes 

V c (t) < 1 + J (1 + T)V c (r) 10 e 4DT exp (v c (r) 9 (l + r) 2 e 4Dr )dr. (4.10) 

By Gronwall's inequality, there exists an interval [0,T) independent of c where V c (t) 
remains finite for all c > 1, i.e., V c (t) < 1- Using this estimate we can complete the 
proof of Theorem [3J Let T^ ax denote the maximal time of existence of a solution of 
H1.3fl - 11.4f> and assume T^ lax < T for some c > 1. Since V c (t) is an increasing function 
of time, this implies V c {T^ a ^) < oo and so, by Theorem^ T^ ax — oo, a contradiction. 
Hence, the solution is defined on the interval [0,T) for all c > 1. 

Proof of Theorem For brevity we omit stating that the estimates below are valid 
for t € [0,T']. The claim on d t 4> follows directly by (|4.7|l and the assumption V c {t) < 
1. Note also that, by l|4.5[) . V x 4> — C( c ~ 2 )- We start with estimating <j>. By the 
assumption we made on <j?± and Lemma |21 the second term in (|3.2|) becomes 



t 

47T 



<pY(x + etui) du = C(c" 3 ). (4.11) 
M=i 



For the first term in (|3.2|) . with the assumption on <fi Q n , Lemma |3 and the fact that 
Ag$ = 4?r J / in dp, we get 

9 *f^ / $(x + ctLj)du) =-(4 7 r)- 1 c- 2 / ^4dy 

= "W f^dpdy. (4.12) 

J\y-x\>ctJ \y-X\ 
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For i/)(t,x) in (|3.1[) . first consider 

f(t - c~% - x\,y,p) f(t - c~ x \y - x\,y,p) 



y — x I < < ' 



yl-j-c 2 p 2 \y — x\ 



dpdy 

< -/ ; c~ 2 f(t- c^\y- x\,y,p)dp dy 

\v\<R+r c (t)t J \p\<v c (t) \V ~ x \ 



<rr\ 



Then ?/> becomes 

V'(i) a:) = — cT 2 
So collecting l|4.11|l . (j4.12fl and l|4.13[) . we get 



f(t - c-'ly - x\,y,p) dpr-^-r + O^ 4 ). 

|y-x|<ctJ I?/ — ^1 



(4.13) 



y — x \<ct 



f(t~c 1 \y - x\,y,p)dp 



dy 



\y-x\ 



\y-x\>ct 



\y - x \ 



/(max{0,t-c x |y- x|},j/,p)dp 



rfiy 



+ C( C - 3 ). (4.14) 



So using l|2.8(l and Q4.14fl . now we estimate 



\c 2 (j)(t,x) - U(t,x)\ 



0(c- 1 )+ // f 00 (t,y,p)dp 



dy 



\y-x\ 



/(max{0,t-c 1 \y - x\},y,p)dp- 



dy 



< 



\y-x\ 

/oo(max{0,t- c~ x \y- x\},y,p) - foo(t,y,p) 



dp 



dy 



\y-x\ 



/(max{0, t-c l \y- x\},y,p) 
/oo(max{0,f- c~ l \y- x\},y,p) 



a d y , -i 

dp- r + c 



\y-A 



(4.15) 



Define 



Dir(t) :=sup{|/(T,a;,p)-/ 00 (T,x,p)| :t<e [0,i],ieR 3 andpeR 3 }. 
Also define V x (t) as the following 

"Poo(i) = sup {|p| : (z,p) <E supp /oo(s)} + 1. 

0<s<t 

Since (foo, U) is a C 1 solution of (|2.7() - (|2.9|l and the initial data /oo has compact 
support, T^oo is well defined for all t > 0. Note that <9t/oo is bounded on R 6 x [0,T']. 
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Also let V(t) := V c {t) +7 3 00 (t). Then (|¥TT31> becomes 
|c 2 0(t,a;) - U(t,x)\ 
< 



\dtfoo(s, y,p)\ dsdp- ° J 



y|<-R+'P(t)t ^ H^T^t) Jmax{0,t-c- 1 |j/-x|} 12/ — x \ 



dy 



[ [ DpimaxiO^-c-^y-x^dp- 

J\y\<R+V(t)tJ\p\<V(t) W - x \ 



c-\ 



^Dp^+c- 1 . (4.16) 

Before estimating Dp(t), let us look at V ' x (f>. In the process of proving Theorem [3 
we have also shown that 

d x ^ om - c-H-i S ]y _ x]=ct J (1+c ^ 1+c ^ (v,p) d P dS y + 4, 

= -c- 2 // ^/(max{0, t - c- 1 !?, - + 0(c^), (4.17) 

cf. JOJ . For II Xi in we use S(</>) = ©(c^ 1 ) and < 7> c (t) 4 . Therefore 

\H Xi \<c- 3 [ f c- 1 f(t-c- 1 \y~x\,y,p)dp-^L_ <c-\ (4.18) 

j|i/i<i?,+-p c (t)t "'ipi<'Po(t) ly x i 

The estimation of //X^ is similar to the one of II Xi - Recall that |c"^ x i | < V c (t) 4 
and V x 4> = 0(c~ 2 ). So we get 

\IH x A<c- 4 . (4.19) 
Now collecting Ij4.17|l - H4.19j) we obtain 

d Xi 4> = -c- 2 jj{Vi- Zi)/(max{0, i - c" 1 ^ - x \}, y,p) dp r^r 3 + O^). (4.20) 

By the similar argument in 14.16fl . with l|2.8|) and 14.20fl . we estimate 
\c 2 V x 0(t,x)-V x U(t,x)\ 

0(c- 1 )+ f[(y-x)f 00 (t,y,p)dp " !! 



\y - x\ 3 

ff( y ~ - T )/( max {0) t - crl \y - %\},y,p) d P \ y d ^ x v. 



/ / / \dtfoo(s,y,p)\ dsdp 

J\y\<R+V(t)t J\p\<V(i) Jma.x{0,t~c- 1 \y-x\} 

I f DFim&xiO^-c-^y-xWdp— 

J\v\<R+V(t)tJ\v\<V(t) \V 



\y - x\ : 



<D F (t)+c- L . (4.21) 
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To estimate Dp(t), let us define Df := f — f^. Then using the two Vlasov equations 
l|1.4ll and (|2.7|l . we obtain 



d t D f +p-V x D f - [S(<f>)p 

= (p-p)-V B /oo+ [S(0P 



c 2 V x <; 



y/1 + C~ 2 p 2 

c 2 V x d 



y/l + C~ 2 p 



] ■ V P £>/ (4.22) 
V K C7] •V p / 00 +45(^)(/ 00 + £>/). 



Note that \p — p\ < c 2, P c (t) 3 . Also note that V^/oo, V p /oo and foo are bounded on 



[0,T']. Then with JOTJ, (|4~2ll) becomes 



d t D f +p-V x D f - [S(4>)p- 



sj\ + C- V 



] ■ V p Z>j 



<r.- x 



c 2 - V x [/ 1 + c- 1 1 £>, 1 < (i ) 



(4.23) 



Using the characteristics (X,P)(s) of i|1.4(l with (X,P)(t) = (ir,p), compute 

^(s.XW,^))! <D F (,s)+c- 1 . 
Note that £>/(0, AT(0), P(0)) = 0. Therefore integrating lET2lj) we get 



D F(t) < [ D F (s)ds + c- 1 . 
Jo 



(4.24) 



So GronwalPs inequality implies Dp(t) < c 1 , which gives l|2.14|l . The proof of (|2.13|) 
is completed by (I4.1(i|) and 14.2111 . This concludes the proof of Theorem [3] 
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